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We present the latest results from the UKQCD/RBC collaborations for the A73 form factor with 
2 + 1 flavours of dynamical domain wall quarks. Simulations are performed on 16 3 x 32 x 16 and 



24 3 x 64 x 16 lattices with three values of the light quark mass, allowing for an extrapolation to 
the chiral limit. After interpolating to zero momentum transfer, we obtain the preliminary result 
f+ n (0) — 0.9680(16), which is in excellent agreement with an earlier Nf = 2 result. 
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1. Introduction 

— > 7lZv (A73) decays provide an excellent avenue for an accurate determination of the 
Cabibbo-Kobayashi-Maskawa (CKM) ^ quark mixing matrix element, \V US \. This is done by 
observing that the decay amplitude is proportional to |Ki| 2 |/+(<? 2 )| 2 > where f+(q 2 ) is the form 
factor defined from the K —> % matrix element of the weak vector current, Vp = sy^u 

(n{p')\V„\Kip)) = (p^+P^)f + (q 2 ) + (p^-p^f-^ 2 ), q 2 = (P-P'f • (1.1) 

In chiral perturbation theory (ChPT), /+(0) is expanded in terms of the light pseudoscalar 
meson masses, m K , mx, 

/+(0)=l+/ 2 +/4 + ..., (/» = <?«,*,„))■ d-2) 

Conservation of the vector current ensures that /+(0) = 1 in the SU (3) flavour limit, while SU (3) 
flavour breaking effects occur only at second order in (m s — m U( f) due to the Ademollo-Gatto Theo- 
rem [Q], which states that f% receives no contribution from local operators appearing in the effective 
theory. As a result, fi can be determined unambiguously in terms of m n , and f n , and takes the 
value f2 = —0.023 at the physical masses 
Our task is now reduced to one of finding 

A/ = /+(0)-(l + / 2 ). (1.3) 

In order to obtain a result for / + (0) which is accurate to ~ 1%, it is sufficient to have a 20-30% error 
on Af. Until recently, the standard estimate of Af = —0.016(8) was due to Leutwyler & Roos [0], 
however a more recent ChPT analysis favours a positive value, Af = 0.007(12) A calculation 
of Af on the lattice is therefore essential. 

The last few years have seen an improvement in the accuracy of lattice calculations of this 
quantity [|5| §,0, §], with the results favouring a negative value for Af in agreement with Leutwyler 
& Roos. The most recent study used 2 flavours of dynamical domain wall fermions to obtain a 
result Af = -0.009(9) [§]. 

The UKQCD and RBC collaborations have embarked on a program to improve on earlier 
studies by using Nf = 2 + 1 flavours of dynamical domain wall fermions at light quark masses and 
on large volumes. We present here preliminary results from this study. 



2. Lattice Techniques 
2.1 Parameters 



We simulate with Nf = 2+1 dynamical flavours generated with the Iwasaki gauge action [J10J] 
at p = 2.13, which corresponds to an inverse lattice spacing a -1 1.6GeV [11], and the domain 
wall fermion action [|l2|] with domain wall height M$ = 1.8 and fifth dimension length L s = 16. 
This results in a residual mass of cu71xq$ ^ 0.00308(3) [ ]1 lp . The simulated strange quark mass, 
am s = 0.04, is very close to it's physical value JTT|], and we choose three values for the light quark 
masses, am u d = 0.03, 0.02, 0.01, which correspond to pion masses m % f» 630, 520, 390 MeV [p"T|]. 
The calculations are performed on two volumes, 16 3 x 32 and 24 3 x 64, at each quark mass. For 
more simulation details, see [JTT1] . 
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Figure 1: Ratio for /o(<7max)' ^ (*'>0> as defined in Eq. (2.2), for three simulated light masses am 



'hid 



0.03, 0.02, 0.01 for two different volumes, 16 x 32 (left) and 24 x 64 (right). Further simulation parameters 
can be found in JTT|]. 

2 - 2 /o(<7max) 

We start by rewriting the vector form factors given in Eq. dOb to define the scalar form factor 

,2 



/o(<f)=/+(<z 2 ) + 



r 



/-(<T 



(2.1) 



which can be obtained at = (m K — m n ) with high precision from the following ratio [13] 



, Cf» (;V;0,0)Cf (t',t;6,6) 



4m K m n 



Mlmax)\ 



cf (fV;0,o)cf(«',/;0,o) i.e-o- 

where the three-point function is defined as 

cJe(^^^,/ ! p)=^e- ^ ^-^e-^<0|^|Q><Q|y M |p)<p|4|0 



(2.2) 



(2.3) 



.v,i' 



with P,Q = n or K and ^W/f) is an interpolating operator for a pion(kaon). We note that R l (t',t) = 1 
in the SU (3) flavour symmetric limit, hence any deviations from unity are purely due to SU (3)fl av0 ur 
symmetry breaking effects. 

In the left (right) plot of Fig. |j] we display our results for R (t' ,t) for each of the simulated 
quark masses as obtained on the 16 3 x 32 (24 3 x 64) lattices. It is immediately obvious that R l (t',t) 
can be measured with a very high level of statistical accuracy. We also note that the ratio becomes 
larger the further we move away from the SU (3)fl av0 ur limit. Since there is no spatial momentum 
involved in this ratio, the results obtained on the two different volumes should agree, and any 
difference can only be due to finite size effects. The two plots in Fig. [l] indicate that within statistical 
errors, finite size effects on /o(<7m ax ) we negligible. 

Finally, we note that the increased time extent of the 24 3 x 64 lattice allows for a longer plateau 
from which we can extract our result. 

2.3 Investigating the momentum transfer dependence 

To study the q 2 dependence of fo(q 2 ), we construct the second ratio 



R 2 (t',t;p',p) 



Cf n {t' ,t;p' ,p)C K {t;0)C n {t' -t-0) 

Cf K (t' ,t;d,0)C K (f,p)C n (t' - t;p') f,(f-t)- 



E K {p)+E n {p') 



F(p',p), (2.4) 
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Figure 2: Ratio for F[p,p'), R (t',t; p',p), as defined in Eq. (2.4), for bare quark mass am uc [ = 0.02 with 



momentum transfer, |ag| 2 = 1 for two different volumes, 16 3 x 32 (left) and 24 3 x 64 (right), 
from which we are able to extract f+(q ) via 



F(jp',p) 



Mi 2 ) 

/o(« 



2 ) 

max/ 



E K (p)-E„(p') 2 , 
1 + ^+^') SW ' 



/+(<7 2 



(2.5) 



in Eq. ( |2.4| ) is the standard pion (kaon) two point function. 

Figure [2| displays a typical example of R 2 (t' ,t;p' ,p) for bare light quark mass am u< i = 0.02 and 
momentum transfer \aq\ 2 = 1, where we have set p' = and averaged over equivalent 3-momenta 
corresponding to the same 4-momentum transfer. The left plot shows our result from the 16 3 x 32 
lattice, while the right displays 24 3 x 64. Note that, unlike in Section 2.2, we are now including 
finite spatial momentum, so disagreement between the results on the two different volumes is not 
an indication of finite volume effects in this case. 

Before we can extract f+{q 2 ) from Eq. (2.5), we need to calculate ^(q 2 ) = f-{q 2 )/ f+(q 2 )- 
This is achieved by constructing a third ratio: 



R 3 k (t',t;p',p) 



We then obtain ^ (q 2 ) from 



C^(t',t;p',p)Cf K (t',t;p',p) 
Cf*{t',t;p>,p)Cf K {ti,t;pi,p) 



(£=1,2,3) . 



(2.6) 



-(E K {p) + E K {p')){p + p') k + (E K {p) +E K (p')){p + p') k R\ 
(E K (p) + E K (p>) ) (p - p% - (E K (p) -E n {p')){p + p') k R\ 



(2.7) 



We observe that R^.(t',t;p',p) = 1 in the S£/(3)fl avoU r symmetric limit, and deviates only 
slightly from unity at our simulation quark masses. Consequently, from Eq. ( |2.7[ ), £,{q 2 ) has a 
small magnitude ^0.1 with an error typically 25%-100%. 

Finally, we can double the number of available q 2 values by repeating the steps above (Eq. (2.4)- 
(p?7|)) for the % — > K matrix element as described in Section V of Ref. []|]. 

3. Results 

3.1 Interpolation to q 2 = 

We are now in a position to combine the results obtained above for the /o(g ma x)> F(p,p') and 
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(aqr 



Figure 3: Scalar form factor fo(q 2 ) for bare quark masses am uc i — 0.03 (left) and 0.02 (right). Results are 
obtained on two volumes V = 16 3 x 32 (red diamonds) and V = 24 3 x 64 (blue triangles). The solid line in 



each plot is the result of a fit using a monopole ansatz (Eq. (3.2)) 



E, (q 2 ) to reconstruct the scalar form factor 



nn 



(3.1) 



We present our results obtained on each volumes for fo(q 2 ) in Fig. ||for quark masses am u d = 0.03 
(left) and am U( j = 0.02 (right). In the intermediate q 2 -range in both plots, we see good agreement 
between the results obtained on the two different volumes, indicating that finite size effects are 
negligible, at least for the quark masses considered here. This means that we now have results over 
a large range of q 2 to fit to. 

We fit our data with a monopole ansatz 



Hq 2 



MO) 



(l-q 2 /M 2 ) ' 



(3.2) 



which we find to describe our data very well. This enables us to interpolate our results for fo(q 2 ) 
and /o(<7max) to q 2 = 0- Future work will involve testing several ansatze in order to obtain /o(0), 
although previous work suggests that results at the quark masses considered here are insensitive to 
the choice of interpolating method. 

3.2 Chiral Extrapolation 

Now that we have obtained results for /+(0) = /n(0) at three different quark masses, we are 
in a position to attempt an extrapolation to the physical pion mass. Inserting these results into the 
expression given in Eq. (|T"3|), together with f2 calculated at the simulated quark masses using the 
ChPT formula @, [l|], we are now left with the task of chirally extrapolating Af. 

The Ademollo-Gatto Theorem implies that Af oc (m s — m ut f) 2 , hence we attempt a chiral ex- 
trapolation using 



Af = a +B(m s — m i 



ud) 



(3.3) 



Note that in the SU (3)fl avour limit, Af = 0, so we expect that a fit to our data should produce a « 0. 
In the left plot of Fig. H, we show the chiral extrapolation of our results using a slightly modi- 



fied version of Eq. (3.3) which allows the result in the chiral limit to be obtained from the intercept 
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Figure 4: The left plot shows the chiral extrapolation of A/ using a trivial modification of Eq. (3.3). The ver- 
tical dotted line indicates the SU (3)fl av0 ur limit. The right plot is an alternative chiral extrapolation (Eq. (3.4)) 



with the y-axis. We are encouraged by the fact that Af passes through zero at the SU (3)fl avour sym- 
metric point (denoted by the vertical dotted line), and we find in the chiral limit Af = — 0.0090(1 1). 
Alternatively, it has been noted that is convenient to consider an extrapolation of the ratio | 



R 



A/ 



rat 



-m 2 ) 2 



a + b{m K + m n ) 



(3.4) 



Extrapolating our data using this form provides an estimate of the systematic error in the choice 



of chiral extrapolation ( |3.3| ). This obviously requires further investigation and will be reported on 
in a forthcoming publication. We show the extrapolation using Eq. (3^) in the right plot of Fig |] 
from which we extract a result at the physical meson masses (vertical dotted line). We take the 
difference between the results obtained from the two extrapolations (0.0011) as an estimate of the 
systematic error due to the chiral extrapolation. 

Finally, since we only have results at one lattice spacing, we are unable to extrapolate to 
the continuum limit. However, lattice artefacts are formally of <^(a 2 AQ CD ) w 2.5%. Hence our 
preliminary result is 



Af = —0.0090(11) (11) (2) 



fl n {0) 



: 0.9680(16) , 



(3.5) 



where the first error is statistical, and the second and third are estimates of the systematic errors 
due to the chiral extrapolation and lattice aritefacts, respectively. This result agrees very well with 
a recent two-flavour result [^], also obtained with domain wall fermions at a similar lattice spacing 
(a~ l sa 1.64 GeV), indicating that the effects due to a dynamical strange quark are small. 
Using |V U i/ + (0)| = 0.2169(9) from the experimental decay amplitude [15]: 



we ftnd\V ud \ 2 + \V us \ 2 + \V llh \ 2 
by PDG(2006) 8 = 0.0008(1 V 



\V US | = 0.2241 (9) exp (4) /+(0) (3.6) 
1 — 5, 8 = 0.0015(7) which can be compared with result given 



4. Summary and future work 

We have presented a preliminary result for Af = /+(0) — (1 + fi) using Nf • = 2 + 1 dy- 
namical domain wall fermions with three choices for the light quark masses. Our result Af = 
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—0.0090(11) (11) (2) agrees very well with the Nf = 2 result [Q] and confirms the trend of other 
lattice results ^, 0, which prefer a negative value for Af, in agreement with the early result of 
Leutwyler & Roos [|[]. We performed our simulations with matched parameters on two volumes 
and we observe no obvious finite size effects. 

This result can be improved by decreasing the error on the point at am U( i = 0.01 and simulating 
at lighter quark masses. Additionally, this result has been obtained at a single value of the lattice 
spacing, so future simulations will need to be performed at least at one more lattice spacing to 
investigate scaling behaviour. 
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